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' Abstract 

' We discuss time evolution of some solitary waves described in the first part of this work. The 

adiabatic motion of the non-linear non-dispersive waves composed of three lumps is interpreted as 

C ' three-body low energy scattering of these particle-like kinks. 

^3 



; 1 Introduction 

After a brief summary of the stability properties of the solitary waves discovered in the first part we shall 
I devote this second part to the analysis of low energy dynamics. Study of Manton's adiabatic motion of 
' solitary waves [22\ is very simple in this case because the geometric metric of the moduli (parameter) 
space that governs the dynamics does not depend on the translational parameter xq; thus, it is always 
possible to apply a transformation that leads to an Euclidean metric. On the other hand, there are only a 
few works that have addressed research into solitary waves in three-component scalar field theory models, 
' see \'2^\ I24| I25L I26j . In this richer case there are three-parametric families of kinks, and the adiabatic 
Q^ ! evolution of three-body lumps is associated with metrics with curvature. The low-energy dynamics of, 

I in this case, three-body solitary waves is therefore much more intricate. 
^Sj ' In this paper, we shall address a three-component scalar field model that generalizes the two- 

■ component model discussed in Reference |14j to three fields. The organization of the paper is as follows: 
, In Section 1 some remarks about the stability of these solutions are given. Finally, in Section 2 we shall 
describe the adiabatic evolution on the moduli space of solutions of configurations composed of three 
basic lumps. 

; 2 Some remarks about the stability of solitary waves 

X 

^ . In order to analyze the stability of the solitary wave solutions, we shall apply different procedures. The 
direct approach is study of the spectrum of the second-order fluctuation (Hessian) operator -formula (10) 
in part I-, now around kink solutions. A non-negative spectrum of the kink Hessian operator ensures 
that a particular solitary wave solution is stable. Zero eigenvalues in the spectrum, zero modes, show 
that certain perturbations send static solutions into static solutions. Thus, zero modes are associated 
to neutral equilibrium directions in a sub-space of the configuration space where the adiabatic evolution 
of solitary waves takes place. It is easy to check that the existence of a kink family 4){x, c) bears the 
presence of an eigenfunction ^ that belongs to the kernel of H, i.e., ^ is a zero mode. When detailed 
information about the spectrum of the Hessian is lacking, an interesting procedure to study the stability 
of kink families is the application of Morse index theorem. A point -referred to as a conjugate or focal 
point- that is crossed by all the members of this family implies the existence of a negative eigenvalues in 
the spectrum of the Hessian operator [HH rTK]. 
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2.1 Stability of basic solitary waves 

Let us consider first basic solitary waves, i.e., the K^^, K^^ and Kf^ kinks. 

1. K'y^ kink. Because only the third component is non-null, the Hessian is a diagonal 3 x 3-matrix 
differential operator with entries: 




The spectral problem is fully solvable because we are dealing with three ordinary Schrodinger 
operators of the Posch-Teller type. Hu and H22 govern orthogonal fluctuations to the kink orbit 

in internal space, whereas "^33 takes into account fluctuations tangent to the orbit. There are 
no discrete eigenvalues in the spectrum of Tin and 7^22 j all the eigenfunctions belonging to the 

continuous spectrum. In the first case, the continuous spectrum starts at the threshold ?|; w^(g') = 

2 2 

Q'^ + ^ is a non-degenerate eigenvalue in the ^] interval but it is doubly degenerate for higher 

3 ^ ^ 4 

values of uP'. In the second case, the threshold is ^ and the doubly degenerate spectrum starts 

at On the other hand, is a discrete eigenfunction of 7^33 with = eigenvalue (zero 
mode) due to spontaneous breaking of translation symmetry by the kink solution. The continuous 
spectrum u'^ = + 2cj| is non-degenerate for uP' G [2c7|, 8c7|] and doubly degenerate for oj^ > 8fT|. 
In conclusion, the lack of non-ncgativc eigenvalues allows us to claim the stability of this solution. 

The Hessian for any other solitary wave solution is a non-diagonal matrix operator and the task of 
solving the spectral problem becomes hopeless. There is, however, a loophole to avoid this problem 
(not of necessary use in the K^^ case). On the (^3 axis, the Hamilton characteristic function 
-formula (25) in part I- becomes: 



The ^-kink orbit is a flow line of the gradient of W , a polynomial function. Thus, the energy is 
given by |VF(0'^) — W{^^)\, a topological bound, the absolute minimum of E in the C^^ topological 
sector. Therefore, K^^ is stable. 

2. i^p^ kink. In this case the Hessian is not diagonal but the alternative argument works exactly in 
the same manner, -formula (25) in part I- now reduces to: 



Again, |P^((?!)'^) — W{(j)^)\ is the absolute minimum of E in the topological sector and the K^^ 
kink is stable. 



W(</.3) = -=[(</.i(x)-l)2 + 2c7i0i(x 



')] 



W{4>iA2) = -)= [(</)! (x) + 4>l{x) - If + 2al^l{x: 



3. K^^ kink. The same argument also works with 

Wi4>2,4>3) = ^ [(M^) + 4>lix) - If + 2al4ix) + 2al4] 
and E{K^^) = \W{4)'^) - M^(0^)|, guaranteeing the stability of K^'^. 
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2.2 Stability of two-body solitary waves 



The general Hessian operator for K2 {b) and K2 {b) kinks is a non-diagonal 3x3 matrix differential 
operator with an unknown spectrum. It is possible, however, to apply the Morse index theorem, see 
[SSISni- Note in Figures 5(c) and 6(c) that there are no conjugate points -points where all the trajectories 
meet- in the congruence of kink trajectories of this type; the Morse index theorem relating the number 
of conjugate points to the number of negative eigenvalues of the Hessian reveals that both the {b) 
and K2^{b) families are formed by stable kinks. The absence of conjugate points is due to the fact that 
the reduced Hamilton characteristic functions -formula (25) in part I- for both kinds of solitary waves 
are regular functions with no conical singularities where the gradient flow of W would be undefined. 

This result is confirmed directly for a particular member, 7^2**^(0), of the first family. In this case the 
Hessian operator -formula (10) in part I- is diagonal, with the following entries: 



and the analysis of the spectral problem is as follows: there are no discrete eigenvalues in the spectrum 

2 2 

of Till. The continuous spectrum uP'{q) = + ^ starts at the threshold ^ and is non-degenerate in the 
[?|, ^] interval but doubly degenerate for higher values. 2^ is a zero mode of "^22 due to spontaneous 

(Jj (Tj OX 

symmetry breaking of translational invariance by the solitary wave. The continuous spectrum = 
+ 2a 2 is non-degenerate for uP' S [2cj2,8cj2] and doubly degenerate for u;^ > Sal. Finally, the third 
component, "Hss, has a zero mode with eigenfunction ^(0) due to motion without energy cost in the 
space of kink orbits. We find a continuous spectrum, which is non degenerate in the range ^1 and 
doubly degenerate from this value. Again all the eigenvalues are non-negative and the stability of this 
solution is established. 

2.3 Stability of three-body solitary waves 

The K^^{a,b) kinks, encompassing the three basic kinks with no static forces between them, form a 
set of distinguished solitary waves of this system. One- and two-body solitary waves arise in this model 
because of the embedding of lower dimensional systems of the same type. They are generic in the sense 
that this variety of solutions depends on as many parameters -two apart from the center of the kinks- 
as possible for a three-component scalar field theory. Moreover, the stability of this two-parametric kink 
family is easily established by a quick and visual application of the Morse index theorem. Away from the 
vacuum points, (^'^ and (f)^, the K^^{a, b) kink orbits do not intersect; there are no conjugate points, see 
Figure 7(c). Being the gradient flow lines of W -formula (25) in part I-, conjugate points cannot exist 
because TV is a regular function and the flow is well defined. Therefore, these solutions are stable. 

Although the Hessian operator is a very complicated 3x3 non-diagonal matrix differential operator 
and although the spectrum is essentially unknown, it is easy to prove the existence of three zero modes: 
1) ^■(x;a,6), 2) ^-(rc;a,6), and 3) ^■(x;a,6). Given a one-parametric family of solutions (j)a{x;a), 
simply differentiate both sides of the static field equations with respect to the parameter: 



W33 



Till 



'H22 




•a 



(x;q) 




dV 
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On the other hand, the spectral problem is solvable for a particular member, the (0, 0) kink, 
living in the (pi axis. The operator is diagonal with entries 

15 

3 tanh V^y — — sech^\/2?/ 
2cr| + 2(t| + (1 - 2cr|) tanh V2y - - sedi'^V2y 
2cr| + 2(t| + (1 - 2cr|) tanh V2y - - sech^ V2y 

where y = log governs the behavior of the fluctuations tangent to the orbit on the solution 

whereas 7^22 and 7^33 govern the orthogonal perturbations to the kink orbit, pushing the solution towards 
another with non-null values of a and b . The discrete spectrum of Tin comprises only a zero mode whose 
eigenfunction is whereas the continuous spectrum u;^ = (7^ + 1 is non-degenerate for o;^ G and 
doubly degenerate for values greater than 4. Likewise 7^22 presents a zero mode with eigenfunction 
1^ and the continuous spectrum is non degenerate in the range [2cr|,2a|]; doubly degeneracy of the 
spectrum starts from this value. Finally, 7^33 holds a zero mode with eigenfunction ^ and [2a^, 2a^] and 
(2 max{cj|, 00) are respectively the ranges of the singly and doubly degenerate continuous spectrum. 
We conclude the stability of this solution. 



7^11 

7^22 

W33 



dy'^ 
dy"^ 
dy"^ 



+ 5 



+ 1 



+ 1 



2.4 Instability of four-body, six-body and seven-body solitary waves 

Note that the K^^{b) and K^'-^{b) kink trajectories always meet at a common point: the conjugate point, 
see Figure 11(c) and 12(c). For the K^^ family, this point is the focus Fi, whereas for the family K^'-' 
it is one of umbilicus points A of E. This is the signal of the instability of these solutions: according to 
the Morse index theorem, the Hessian for any member of these families has a negative eigenvalue |341 135j . 
The undefined flow at the foci of the ellipse ei or the umbilical points of the ellipsoid E arise because 
the sign combinations in the Hamilton characteristic function providing these trajectories are such that, 
back in Cartesian coordinates, is a function which is non-differentiable at either Fi or A. Thus, these 
are flow lines obtained by gluing two flow lines generated by different Hamilton characteristic functions 
at these points . The energy depends not only on the value of the Hamilton function at the beginning 
and the end of the orbit but also on the values at the conjugate points. The absolute minimum value of 
the energy in the corresponding topological sector is non-saturated and four-body kinks are unstable. 

For the one-parametric six-body K^^{b) and K^^{b) kink families, the situation is exactly the same. 
There exists a conjugate point, which is crossed by all the members of each family. For the Kq'~" {b) it is 
the focus F2 of the ellipse 62, and for the KQ^{b) it is the focus F3 of the ellipse 63. Again, application 
of the Morse index theorem and considerations on the impossibility of reaching the absolute minimum of 
the energy when there are conjugate points reveal the instability of six-body lumps. 

The case of seven-body solitary waves ^C^*-" (72, 73) is more complicated. There exist focal lines, curves 
of conjugate points; on each point of the characteristic hyperbola h a whole one-parametric subfamily of 
trajectories meet. The same is true for a point in the characteristic ellipse 64. The Morse index theorem 
tells us that the Hessian for these solutions has two negative eigenvalues. Furthermore, the absolute 
minimum of the energy is not reached by i^^^*-' (72, 73) kinks because of the existence of lines of conjugate 
points crossed by the kink orbits: the instability of seven-body solitary waves is proved. 

The dynamical reason for the instability of these solitary waves is the peculiar configuration of the 
lumps arranged in these solutions. All these solitary waves are static solutions formed by several basic 
particles, where some of these lumps travel superposed. Because they are static solutions, the lumps are 
distributed in order to achieve an equilibrium configuration where no forces between lumps are involved: 
the exact overlapping between lumps produces this balance of forces. Any small perturbation slightly 
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splitting the superposed lumps breaks the unstable equilibrium releasing forces between the lumps. The 
negative eigenvalues of the Hessian operator associated with this process indicate that the evolution of 
these solitary waves moves them away from the initial static configurations. 



3 Three-body low-energy scattering of non-linear waves 

The three parameters xq, a and b of the family of solitary waves fix the center of mass and 

the relative positions of the basic lumps in the composite kink configuration. We shall now study the 
evolution of composite solitary waves within the framework of Manton's adiabatic principle, see [21\ I22j. 
by looking at changes in a and b in time. The Manton adiabatic hypothesis can be summarized as 
follows: the dependence on time of solitary waves develops only in the parameters of the solution, i.e. 
(j)^{x;t) = (l?^{x;a{t),b{t)). This is a good approximation if the kink evolution is so slow that the 
differential equations -formula (5) in part I- are satisfied at every given t with good accuracy. 

The most general family of trajectories arising in a iV-dimensional Hamiltonian system depends on 
2N parameters: separation constants and integration constants. Finite action fixes the separation 
constants at a fixed value. Therefore, the most general family of solitary waves in A^-component scalar 
field theory depends on A^ parameters, the integration constants of the finite action trajectories in the 
analogous mechanical system. In this case adiabatic evolution is of the general form: 

4>^{x;t) = ^'^{x;a\t)) , a = l,2,---,N , I = 1,2,---,N 

Introducing this assumption in the action functional, (??) one obtains the action for a geodesic problem: 



^ — 1 T J 1 



„ . „ „ _ , dt dt 

a=l '' I,J=1 



where the metric tensor gij{a^), the Christoffel symbols ^fj, and the curvature tensor are: 

dx- 



n „ 

RlJKM = ^ / 
i=l ~ 



oo 

dx 



a=l-'-'^ K=l i=l 



RlJKM = —RjIKM = —RlJMK = RjIMK = RrMIJ 



RlJKM + RiKMJ + RiMJK = 



Study of the evolution of solitary waves thus turns into a geometric problem in the adiabatic regime: 
the analysis of geodesic motion in the moduli space of parameters where a non-Euclidean metric is 
inherited from the dynamics of the zero modes. The geodesies obey the associated Lagrange equations: 

d^a^ da^ da-^ 

i,J=i 

Back in our A^ = 3 scalar field theory we consider and ^ = 6 as local coordinates in 

the moduli space of parameters for K^^{a,b) three-body solitary waves. Both the dynamics and the 
geometry of the center of mass determined by the xq parameter are trivial: the component gxQxo is 
constant, the geodesies being straight lines xq = vt + d. Therefore, the dynamics of the center of mass 
obeys Galilean, rather than Lorentzian, velocity transformations in the adiabatic regime. Choosing the 
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values of the coupling constants (t| = 2(t| = ^ and defining the variable z = e ^ as in Section §3, we 
obtain the following integral expressions for the components of the metric tensor: 

= L 2^2(1 + a2. + 52.^ + .3)3 (2) 

, ^, p -ate3(a2 + 2622 + 3.2) 

= L 2V2(l + a2. + 62.2 + .3)3^^ (3) 

, p .[(a4 + 362)^2 + (l + 23)2 + 2a2z(l + 62.2 + ^3)1 

= io 2V2(l + a2. + 62.^ + .3)3 ■ 

Integrals of rational functions in @, (jSJ and @ would be feasible but, depending on the arbitrary 
parameters a and 6, the quadratures involved in this calculation lead to expressions that are so complicated 
that there is no advantage in addressing the problem analytically. Numerical integration, however, gives 
us interesting qualitative information about the metric tensor and the scalar curvature. In Figure 16, the 
components of the metric tensor as well as the scalar curvature are displayed as a function of a and 6. 




Figure 16: Metric Components and Scalar Curvature: a) gaa{o,,b), b) gab{a,b), c) gbb{a,b), d) R(a,b). 

The picture of the scalar curvature is especially useful for obtaining clues about geodesic motion. 
In regions where the curvature is positive the geodesies tend to approach each other whereas where 
it is negative the geodesies move away. Therefore, a certain interaction between lumps arises in the 
study of the evolution of solitary waves even in the adiabatic approximation. Recall that K^^{x]a,b) 
solitary waves are stable and the zero modes are the lowest energy eigenfunctions of the Hessian operator; 
low-energy fluctuations are dominated mainly by these channels, certifying the efficacy of the adiabatic 
procedure. For this reason, the adiabatic method is not applicable to unstable configurations, such as 
four-body, K^^{h) or K^^ {b), six-body, KQ^'{b) and KQ^{b), and seven-body, Kf^ {a, 6), solitary waves, 
because negative eigenvalues arise in the spectrum of the Hessian for these kinks. 

To gain more precise knowledge about adiabatic motion we numerically solve the geodesic equation 
for some different choices of the initial conditions a(0), 6(0), d(0) and 6(0). The values of the parameters 
a{t) and b{t) at t = specify the initial static configuration K^^{a, 6) and the soft perturbations on the 
kinks due to the zero modes ^ and ^ are determined by the values of the time-derivatives a{t) and b{t) 
at t = 0. 

We first describe the scattering process when a perturbation due to the zero mode ^ is exerted on a 
^3(0,6) kink with |6| large. The initial solitary wave starts from some point in the green zone C of the 
moduli space (Figure 9) so that the configuration is formed by two lumps: a K'^^ kink is on the right of 
the energy density profile and the superposition of a Kf'" and a i^f ^ kink is on the left. These lumps 
approach each other, later crash into each other in the red zone A, and finally bounce back, recovering 
their original shape again in the other branch of green zone C, see Figure 17(a). 

The process illustrated in Figure 17(b) is similar to the previous one. Here, a perturbation due to 
the zero mode ^ is applied to a -fC3(a,0) solitary wave. The initial configuration is formed by a K^^ 
lump and the superposition of a K^^ and a K^^ kink. The process of scattering is as follows: the two 
lumps approach each other in the green zone B, collide and travel superposed along the red zone A for 



6 



a period of time and finally split into two lumps, recovering the initial configuration in the other branch 
of green zone B. 




Figure 17: Evolution of solitary waves. 

Although one of the two lumps is a composite of two basic particles, the phenomena described above 
are two-body scattering processes, essentially identical to the collisions of solitary waves studied in j2()j 
and More sophisticated scattering processes are displayed in Figure 18, triggered by perturbations 
due to the two zero modes . Genuine three-body low energy scattering arises, and as a result the main 
novelty in this paper emerges. 

First, start from a K^^{W,-70) solitary wave in the blue zone D near the border with the green 
zone C. In this configuration, the ETp^ and K^'^ kinks are close but the K^^ lump is far away from the 
other two, see Figure 18(a). A perturbation characterized by the values 6(0) = —3, d(0) = 6 is exerted on 
this arrangement of lumps. The evolution of the solitary wave -unveiled from the corresponding geodesic 
found numerically- is as follows: the two closer particles collide and merge whereas the other particle goes 
towards them, all entering together in green zone C. The united lumps travel together and eventually 
split into the initial basic lumps -back again in a blue zone D of moduli space-, the ^ moving towards 
the left and the K^^-^ towards the right on the line. Later, this lump bumps into the K^^ kink, which 
was traveling towards the left. The K^'^ and the K^^ basic particles merge and move as a single one, 
-now in a green zone B-, and, after a certain time, they split, the whole solitary wave again entering in 
other blue zone D in the moduli space. Then, the K^^ kink advances towards the right and the K^'-" 
advances towards the left on the line. Now this particle is attracted by the K^^ lump and they merge, 
coming again to a green zone C. In conclusion, the process can be interpreted as a interchange of the 
particle Kf^ between the other particles K^^ and K^^. 

A last example of three-body collisions is displayed in Figure 18(b). A perturbation determined by 
6(0) = —3, d(0) = 6 is applied on a initial i^3'^(50, — 50) solitary wave living in a blue zone D of the 
moduli space. The K^^{50,-50) confi guration is formed by three separated basic particles. Note the 
apparent attractive forces between the lumps in Figure 18(b) showing the evolution of this solitary wave 
when the perturbation chosen above acts on it. The basic particles tend to approach each other at the 
beginning of the geodesic motion. More precisely, first the K^^ and K^'^ lumps amalgamate when the 
geodesic enters a green zone B of the moduli space, then they travel together and, later on, they all 
merge together with the K^'^ kink. During a period of time, when the geodesic passes through the red 
zone A, the three particles move together and, finally, a process of fission takes place in two stages. The 
K^^ kink splits from the global lump when the geodesic leaves the red zone to enter into the other green 
zone B, and later on so do the K^^ and K^'-" kinks, the geodesic ending back in other blue zone D. 
This process may be repeated if the initial conditions setting the magnitude of the perturbation and the 
initial solitary wave are chosen adequately. 
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Figure 18: Evolution of solitary waves. 
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